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ABSTRACT 

We calculate the nucleon and nuclear photoproduction cross sections for heavy quarks within 
the fc^-factorization formalism, considering the current high energy approaches which include 
nuclear and saturation effects. Our results demonstrate that a future experimental analysis of 
this process would allow to constraint the QCD dynamics at high energies. 



1 Introduction 



The electron-proton [ep) collider at HERA has opened up a new kinematic regime in the study 
of the deep structure of the proton and, in general, of hadronic interactions. This regime is 
characterized by small values of the Bjorken variable x = Q'^/s, where is the momentum 
transfer and ^/s is the center-of-mass energy. The experimental results from HERA have shown 
a striking rise of the proton structure function F2{x, Q"^) for values x < 10~^, which implies that 
the cross section increases faster than logarithmically on energy. This high energy behavior in 
the hard-scattering regime is expected if the underlying dynamics is driven by self-interacting 
massless vector bosons, the gluons. Thus, the steep rise of F2 certainly confirms one of the basic 
predictions of perturbative QCD. As the predictions of the collinear and the /c^-factorization 
approaches, relying on different assumptions, agree with the measurements of the inclusive 
quantities, one current open question is: what is the correct QCD dynamics at high energies? 

In the collinear factorization approach [1] all partons involved are assumed to be on mass 
shell, carrying only longitudinal momenta, and their transverse momenta are neglected in the 
QCD matrix elements. Moreover, the cross sections for the QCD subprocess are usually calcu- 
lated in the leading order (LO), as well as in the next-to-leading order (NLO). In particular, 
the cross sections involving incoming hadrons are given, at all orders, by the convolution of 
intrinsically non-perturbative (but universal) quantities - the parton densities - with perturba- 
tively calculable hard matrix elements, which are process dependent. The conventional gluon 
distribution g{x, fi"^), which drives the behavior of the observables at high energies, corresponds 
to the density of gluons in the proton having a longitudinal momentum fraction x at the fac- 
torization scale /X. This distribution satisfies the DGLAP evolution in /x^ and does not contain 



information about the transverse momenta kj_ of the gluon. On the other hand, in the large 
energy (small-a;) hmit, we have that the characteristic scale fj, of the hard subprocess of parton 
scattering is much less than i/i, but greater than the Aqcd parameter. In this limit, the effects 
of the finite transverse momenta of the incoming partons become important, and the factoriza- 
tion must be generalized, implying that the cross sections arc now A';^-factorizcd into an off-shell 
partonic cross section and a A;_L-unintcgratcd parton density function J-'{x,k±), characterizing 
the /c^-factorization approach [2, 3, 4]. The function is obtained as a solution of the evolution 
equation associated to the dynamics that governs the QCD at high energies. 

Recently, several authors have considered the A;_L-factorization approach in order to analyze 
some non-inclusive observables and they have obtained a better description of these quantities 
than the coUinear approach (For a recent review see Ref. [5]). However, the current situation 
is still not satisfactory, due to the large uncertainty associated to the lack of a complete knowl- 
edge of the unintegrated gluon distribution. A search for this distribution in the nuclcon has 
been subject of active both theoretical and phenomenological research in recent years. One 
of the most promising process to constraint tha quantity is the heavy quark production in jp 
interactions [6, 7, 8]. At high energies, the production of open-flavored QQ pairs is described 
in terms of photon-gluon fusion mechanism, which is expressed in terms of the unintegrated 
gluon distribution and the parton-lcvcl matrix elements. As the behavior of the unintegrated 
gluon distribution is determined by the underlying dynamics, in Ref. [7] the deviations between 
the results obtained using a distribution derived from the saturation model and the DGLAP 
evolution equations were analyzed in detail, as well as compared with the predictions of the 
coUinear approach. Those predictions reasonably describe the current HERA data, with siz- 
able deviations appearing basically at higher energies. These results motivate us to analyze 
the heavy quark production in the context of photonuclear interactions, where the deviations 
among the approaches should be amplified due to the presence of the nuclear medium. Here, we 
estimate the cross section for the nuclear photoproduction of heavy quarks, considering the k_i- 
factorization approach and distinct nuclear unintegrated gluon distributions. In particular, we 
consider the unintegrated gluon distribution obtained recently in Ref. [9], where the saturation 
model was extended to the nuclear case using the Glauber-Gribov formalism. Moreover, we 
also consider the derivative of the EKS nuclear gluon distribution [10], which takes into account 
the nuclear medium effects (shadowing, antishadowing, EMC and Fermi motion). Furthermore, 
we compare our predictions with those obtained from the coUinear factorization approach. In 
this case, we consider the EKS nuclear gluon distribution as input in our calculations. We also 
investigate the possibility of the nuclear gluon distribution to be modified by the high parton 
density effects, as estimated in Ref. [11]. This procedure allow us to estimate the theoretical 
uncertainty present in the predictions of the cross sections at high energies. 

This paper is organized as follows. In the next section we present a brief review of the k±- 
factorization approach for the heavy quark photoproduction and estimate the cross section for 
this process in the nucleon and nuclear cases. In particular, we present the proton unintegrated 
gluon distribution obtained from the saturation model and the derivative of the coUinear gluon 
distribution. The results for the total cross sections are compared with the HERA experimen- 
tal measurements. Similarly, we derive the nuclear unintegrated gluon distribution from the 
extension of the saturation model to the nuclear case and the corresponding quantity using 
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the EKS nuclear gluon distribution. Moreover, we present our predictions for the heavy quark 
nuclear photoproduction cross section in the /c^-factorization approach and compare them with 
those ones from the coUinear factorization approach. Finally, in Section 3 we summarize the 
main conclusions. 

2 Nuclear heavy quark photoproduction in the k± - fac- 
torization approach 

Let us start introducing the theoretical prediction for the heavy quark photoproduction in the 
A;_L-factorization approach for the nucleon (proton) case. The relevant diagrams are considered 
with the virtualities and polarizations of the initial partons, taking into account the transverse 
momenta of the incident partons. The processes are described through the convolution of 
off-shell matrix elements with the unintegrated parton distribution, J^{x, k±) (for a recent 
review, see [5]). The latter can recover the usual parton distributions in the double logarithmic 
limit by its integration over the transverse momentum of the k± exchanged gluon. The gluon 
longitudinal momentum fraction is related to the c.m.s. energy, W^p, in the photoproduction 
case as x = Amq/W^p. A sizeable piece of the NLO and some of the NNLO corrections to the 
LO contributions on the coUinear approach, related to the contribution of non-zero transverse 
momenta of the incident partons, are already included in the LO contribution within the k±- 
factorization approach. Moreover, the coefficient functions and the splitting functions giving 
the coUinear parton distributions are supplemented by all-order asln{l/x) resummation at 
high energies [12]. Consequently, in principle, we expect that in the asymptotic regime of 
large energies the use of the /c^-factorization implies an enhancement of the cross sections in 
comparison with the predictions obtained in the coUinear factorization [2]. This is associated 
to the opening of the k± phase space (away from the coUinear region). 

The cross section for the heavy-quark photoproduction process is expressed as the convolu- 
tion of the unintegrated gluon function with the off-shell matrix elements, where the LO matrix 
elements are well known in the literature [5, 6, 7]. The expression for the photoproduction total 
cross section considering the direct component of the photon can be written as [7] , 

atriw,,) fdz £p,, d^k. a.{^.^):F{x,kl■,^.^) 

.{,.,a-.n(t-^)^'»«(^+i;)]. w 

where Di = p'lj_ + tjiq and D2 = {k± — Pi±)'^ + rriQ. The transverse momenta of the heavy 
quark (antiquark) are denoted by Pij_ and P2± = {k± — Pi±), respectively. The heavy quark 
longitudinal momentum fraction is labeled by z. The scale /i in the strong coupling constant in 
general is taken to be equal to the gluon virtuality, in close connection with the BLM scheme 
[13]. On the other hand, in the leading ln(l/x) approximation, ct^ should take a constant 
value. In our further analysis on heavy quarks we use the prescription fi'^ = k']_ + uiq. In Eq. 
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Figure 1: The charm and bottom photoproduction cross sections as a function of c.m.s. energy 
W-fp in the proton case. In the plots are shown the results for the semihard approach (solid 
lines), the saturation model (dot-dashed lines) and collinear (long dashed lines). 

(1), the unintegrated gluon function was allowed to depend also on the scale fi^, since some 
parametrizations take this scale into account in the computation of that quantity [5]. 

In order to perform a phenomenological analysis within the fc^-factorization approach, in 
the following we use two well known parameterizations for the unintegrated gluon distribu- 
tion. First, one considers the derivative of the collinear gluon parton distribution function, 
^dgiuon(a;, ^1) = d xG{x, k']_) /dlnk^, where the GRV94 LO parameterization [14] was consid- 
ered for technical simplicity (a comparison using other pdfs can be found in [7]). Hereafter, 
we will denote as semihard the results obtained using this procedure. It is important to em- 
phasize that if more recent parameterizations are considered, a less steep growth on energy is 
obtained, but still in agreement with the experimental data [7]. The low transverse momen- 
tum region was computed using the ansatz as in Ref. [7]. A shortcoming with this simple 
parameterization for J-' is that it becomes negative at large x. This can be overcome through 
the introduction of the doubly logarithmic Sudakov form factor Tg{k']_, jj^), which gives the 
survival probability that the parton with transverse momentum remains intact in the evolu- 
tion up to the factorization scale /i^ [15]. So, we can write the unintegrated gluon distribution 
as J?Mgiuon(a^, ^1) = d[Tg{k']_, fi'^) xG{x,k']_)]/d\nk']_. We will denote as semihard the results 
obtained using this procedure. 

Second, one compute the heavy quark cross sections using the phenomenological saturation 
model [16], whose unintegrated gluon distribution is given by, 

kl) = 7^ (t^) -P (-7^) (1 - , (2) 

Att^Os \Qi{x)j V Qsix)/ 

where one has used the parameters from [16], which include the charm quark with mass rric = 1.5 
GeV. The saturation scale Ql{x) = {x^/x)^ GeV^, gives the onset of the saturation phenomenon 
to the process. The last factor in the equation above takes into account the low energy threshold 
effects. In Fig. 1 the charm and bottom photoproduction cross sections [17, 18, 19, 20] are 
shown comparing the referred parameterizations for the unintegrated gluon function. As already 
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found in Ref. [7], the saturation model gives a lower bound for the cross sections, having a 
mild increasing on W^p , whereas the derivative of the usual pdf presents a steeper growth on 
energy. The latter disagree with the low energy data, but this can be overcome by introducing 
the Sudakov form factor discussed before. In the plots we also present the results considering 
the coUinear approach for the LO process — > QQ. We have performed fully LO calculations, 
including LO parton densities and a one-loop calculation evaluation of a^. Similar results are 
obtained using NLO calculations, since the NLO corrections for heavy quark photoproduction 
cross section can be expressed in terms of K- factors, which is about 1.6 for bottom and ^/s — 1 
TeV, while the NLO parton densities are smaller by a similar factor than the LO one. Here we 
have used rric = 1.5 GeV, m^, = 4.5 GeV, factorization scale /x^ = s and the GRV94 LO gluon 
distribution. If the the GRV98 LO parameterization is used wc have that the cross section is 
reduced by ~ 20%. In Ref. [21] we have analyzed the dependence of the heavy quark cross 
section in the choices of mass, pdfs and factorization scale. We have that the semihard result 
using .Fdgiuon(3:, k\) is quite similar to the coUinear one at high energies and the deviation with 
the saturation model is smaller for the bottom case (For a more detailed discussion see Ref. 



Having reviewed the notions of the /c^-factorization approach applied to the heavy quark 

photoproduction and computed the relevant cross section at the nucleon level, now we will 
consider an attempt of its extension to the nuclear case. In order to do so, we need a parame- 
terization for the nuclear unintegrated gluon function. Corroborated by the good results using 



where xGa{x,Q^) is the nuclear gluon distribution, which was taken from the EKS parame- 
terization [10]. Consequently, with this procedure we include in our calculations the medium 
effects (shadowing, antishadowing, EMC and Fermi motion effects) estimated by this param- 
eterization. Moreover, we emphasize that this nuclear gluon distribution is solution of the 
DGLAP evolution equations, i. e. it is associated to a linear dynamics which does not consider 
dynamical saturation effects. Similarly to the nucleon case, we denote as semihard the results 
obtained using this procedure. 

In comparison, we calculate the nuclear cross section using the model proposed in Ref. [9] , 
which is an extension of the ep saturation model through Glauber-Gribov formalism. In this 
model the cross section for the heavy quark photoproduction on nuclei targets reads as [9, 22] 



where longitudinal contributions are suppressed and the transverse wave function is known (See 
e.g. Ref. [24]). The nuclear dipole cross section is given by. 



[7]). 




(3) 





a^,p{x, r^ A) = ^d^ft 2 { 1 - exp [-^ AT^(6) alp{x, r^)] | , 



(5) 
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where b is the impact parameter of the center of the dipole relative to the center of the nucleus 
and the integrand gives the total dipole-nucleus cross section for fixed impact parameter. The 
nuclear profile function is labeled by T^(6), which will be obtained from the 3-parameter Fermi 
distribution for the nuclear density [23]. The parameterization for the dipole cross section 
takes the eikonal-like form, o-^jp(.x, t^) = (Jo[1 — exp {—Ql{x) ], with the saturation scale 
previously defined and x — {Q^ + Am'Q)/W; 
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The equation above sums up all the multiple elastic rescattering diagrams of the qq pair 
and is justified for large coherence length, where the transverse separation r of partons in the 
multiparton Fock state of the photon becomes as good a conserved quantity as the angular 
momentum, i. e. the size of the pair r becomes eigenvalue of the scattering matrix. It is 
important to emphasize that for very small values of x, other diagrams beyond the multiple 
Pomeron exchange considered in this formalism should contribute {e.g. Pomeron loops) and a 
more general approach for the high density (saturation) regime must be considered. However, 
we believe that this approach allow us to obtain lower limits of the high density effects in the 
RHIC and LHC kinematic range. Therefore, at first glance, the region of applicability of this 
model should be at small values of x, i.e. large coherence length, and for not too high values 
of virtualities, where the implementation of the DGLAP evolution in the ep saturation model 
should be required. Therefore, the approach is quite suitable for the analysis of heavy quark 
photoproduction in the kinematical ranges of the future eA colliders (eRHIC/TESLA). In Ref. 
[22] we have analyzed in detail the behavior of the dipole-nuclei cross section and estimated 
the nuclear heavy quark photoproduction cross section in the inclusive and diffractive cases. In 
particular, we have predicted large cross sections at eRHIC/TESLA energies. 

The corresponding unintegrated gluon distribution can be recovered from a Bessel-Fourier 
transform to the momentum representation [9] , 



nuc 




which depends on the transverse momentum k± through the scaling variable r = k^/Ql. The 
unintegrated gluon vanishes asymptotically at /c^ — > 0, oo and its maximum can be identified 
with the saturation scale Qsa{x) [9, 25]. 

It has been verified in Ref. [22], where one computes heavy quark photoproduction within 
the saturation approach, that the resummation of high density effects at the proton level is 
less sizeable in the final results at nuclear level. Hence, for heavy quark production, we can 
approximate the dipole nucleon cross section as cr^ipix, r^) ~ cr^Q^x) and then compute 
analytically the unintegrated gluon distribution. Following the procedure in Ref. [26] (See also 
Ref. [27]), we obtain that the nuclear unintegrated gluon distribution may be expressed as 

*i. = ^ (q^) exp (-^) , (7) 
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Figure 2: The charm and bottom nuclear photoproduction cross sections as a function of c.m.s. 
energy W . In the plots are shown the results for the semihard approach (solid lines), the 
saturation model (dot-dashed lines) and collinear using the parameterizations for the nuclear 
gluon distribution: EKS (long-dashed lines) and AG (dashed lines). 



where Q^^(x) = \ATAih) Ql{x) define the nuclear saturation scale. This result is consistent 
with that obtained in Ref. [9] in the first scattering approximation. Such an approximation 
is justified in the heavy quark case, which is dominated by small dipole configurations (large 
transverse momentum fc^ ~ mg) (See discussion in Ref. [22]). In the case of light quarks, the 
process receives sizeable contribution from large dipole configurations, i.e. k\ < Q^, and higher 
order rescatterings are needed for the proton as rendered by the phenomenological saturation 
model. It is important to emphasize that the expression (7) shows clearly the scaling property 
on the variable r = fc^/Q^^. This implies geometric scaling on r in the nuclear heavy quark 
production, which it already been shown in the nucleon case [28]. 

In order to compare our predictions for the nuclear photoproduction cross section obtained 
using the fc^-factorization approach with the collinear one, let us present a brief review of the 
latter. In this case the cross section is given by a convolution between the partonic cross section 
for the subprocess ^yg —>■ QQ and the integrated gluon distribution for the nucleus xGa{x, Q'^)- 
Here we consider the EKS and AG parameterizations for this distribution. As discussed before, 
the EKS parameterization was obtained from a global fit of the nuclear experimental data 
using the DGLAP evolution equations, which is a linear evolution equation which does not 
consider dynamical saturation (high density) effects. In Ref. [11] a procedure to include 
these effects in the nuclear gluon distribution was proposed, resulting in a paramerization 
for this distribution (AG parameterization), which also includes those present in the EKS 
parameterization. The main characteristic of this parameterization is that it predicts a stronger 
reduction of the growth of the gluon distribution at small values of x than the EKS one. In order 
to estimate the sensitivity of nuclear heavy quark photoproduction in these effects, we also use 
this parameterization as input in our calculations using the collinear factorization approach. 

In Fig. 2 are shown the results for the charm and bottom photoproduction cross section as 
a function of energy for lead and calcium. We compare the predictions of the fc^-factorization, 
considering the unintegrated gluon distributions discussed above (SAT-MOD and SEMIHARD 
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in the figure) with resuhs of the colhnear factorization and two parameterizations for the nuclear 
gluon distribution (EKS and AG in the figure). We have that the results obtained using the 
derivative of the EKS gluon distribution and the A;_L-factorization (SEMIHARD in the figure) 
are quite similar with those ones coming from the colhnear approach where nuclear effects 
(EKS parameterization [10]) and high density corrections (AG parameterization [11]) arc taken 
into account. In particular, we have that the predictions using the AG parameterization in 
the colhnear approach are similar to the semihard one, which does not consider high density 
effects. This demonstrate that in this process we cannot distinguish if the modification in 
the behavior of the cross section is associated to high density effects in the colhnear approach 
or a generalization of the factorization without high density effects in the unintegrated gluon 
distribution. On the other hand, if these effects are present and the factorization of the cross 
section is given by the A;_L-factorization, as is the case for the predictions from the saturation 
model, we have that the difference between the cross sections is large, which should allow to 
discriminate between the theoretical approaches. Therefore, the nuclear cross section would 
provide a strong test concerning the robustness of the saturation approach in describing the 
observables. The situation is less clear comparing the semihard approach and the colhnear one. 
One possible interpretation for this result is that the expected enhancement in the semihard 
approach, associated to the resummation of the {ag In ^)"' in the coefficient function [2], is not 
sizeable for inclusive quantities in the kinematic region of the future colliders. Probably, a more 
promising quantity to clarify this issue would be the transverse momentum distribution. In 
this case, the semihard approach seems to be in better agreement with experimental data in 
the pp collisions than the colhnear approach [6] . 



3 Summary and Conclusions 

We calculate the nucleon and nuclear cross sections for heavy quark photoproduction within 
the A;j_-factorization approach. Two simple parameterizations for the unintegrated gluon dis- 
tribution in the proton were considered: the derivative of the colhnear gluon pdf and the 
phenomenological saturation model. The latter underestimate the measured cross section at 
HERA regime and provides a lower limit for the production. On the other hand, the former 
agrees with the available high energy data and produces similar results as the colhnear ap- 
proach. In the nuclear case, we have introduced a simple ansatz for the nuclear unintegrated 
gluon function, namely from the nuclear colhnear gluon pdfs. The Glauber-Gribov extension of 
the saturation model to the nuclear collision is also analyzed, where we found an analytic ex- 
pression for the unintegrated gluon function using a non-saturated (r^ — > 0) proton dipole cross 
section. Such a function shows a clear geometric scaling property on the variable r = k'^/Qlj^. 
The nuclear cross section for lead and calcium in the energy range relevant for the future eA 
colliders eRHIC/THERA were computed, comparing them with the results from the colhnear 
approach which take into account nuclear medium effects and also high parton density phe- 
nomenon. Our results indicate that a future experimental analysis of this process can be useful 
to constraint the QCD dynamics at high energies. 
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